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On the distribution of initial masses of stellar clusters 
inferred from synthesis models 
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Abstract The fundamental properties of stellar clus- 
ters, such as the age or the total initial mass in stars, 
are often inferred from population synthesis models. 
The predicted properties are then used to constrain 
the physical mechanisms involved in the formation of 
such clusters in a variety of environments. Population 
synthesis models cannot, however, be applied blindy to 
such systems. We show that synthesis models cannot 
be used in the usual straightforward way to small-mass 
clusters (say, M < few times 10^ ^o)- The reason is 
that the basic hypothesis underlying population synthe- 
sis (a fixed proportionality between the number of stars 
in the different evolutionary phases) is not fulfilled in 
these clusters due to their small number of stars. This 
incomplete sampling of the stellar mass function results 
in a non-gaussian distribution of the mass-luminosity 
ratio for clusters that share the same evolutionary con- 
ditions (age, metallicity and initial stellar mass distri- 
bution function). We review some tests that can be 
carried out a priori to check whether a given cluster 
can be analysed with the fully-sampled standard pop- 
ulation synthesis models, or, on the contrary, a prob- 
abilistic framework must be used. This leads to a re- 
assessment in the estimation of the low-mass tail in the 
distribution function of initial masses of stellar clusters. 
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1 Introduction 

The study of astrophysical objects is often limited by 
our ability to infer their physical properties such as dis- 
tances, masses or ages, from their observed fluxes. In 
stellar astrophysics, when the distance to an observed 
star is known, the luminosity at different wavelengths 
or spectral energy distribution (SED), Ix, provides con- 
straints on the effective temperature and on its mass- 
luminosity ratio. This ratio depends mainly on the ef- 
fective temperature, gravity and the mass of the star. 
In the case of stellar clusters with known distances, the 
observed luminosity is the sum of the luminosities of 
the stars in the cluster, each one with its own mass- 
luminosity ratio. In a cluster composed by N^, stars, 
this luminosity can be written as 



L. 



cluster 



(1) 



As stated in Eq ((T]), the integrated luminosity of the 
cluster does not provide a great deal of information on 
the stars in the cluster. However, we know that the 
possible luminosities and spectral shapes of individual 
stars are in the range defined by stellar evolution, and 
thus it is possible to group individual stars in represen- 
tative classes or evolutionary phases j of luminosity £j . 
Assuming a total number of iVdass stellar evolutionary 
phases, Eq. [T]can be rewritten as 
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(2) 
(3) 



The problem now becomes the estimation of the nj co- 
efficients in such a way that we can obtain physical 
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properties of the cluster from them. At this stage, it is 
not possible to know the total number of stars in the 
cluster from its integrated light, nor how many stars 
are in a given evolutionary phase. However, we can 
relate the relative number of stars in different evolu- 
tionary phases thanks to stellar evolution: the number 
of stars in a given evolutionary phase is proportional 
to the amount of fuel that can be consumed in that 
phase, and therefore with the lifetime tj of the stellar 
evolut ionary phase. This is the Fuel Consumption The 



orem (Tinslev fc Gunn 



197d iRenzini fc Buzzoiiilll986 



Marigo fc Girardill200lf ) which underlies, implicitely or 
explicitely, any population synthesis method. 

The comparison of different evolutionary phases (say, 
phase i vs phase j) provides the number or population 
ratio Ui/uj in the limit of an infinite number of stars 
in the cluster. Indeed, if the cluster actually has a very 
large number of stars in all the theoretical evolutionary 
phases, we have that 
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These relations hold for the post main sequence evo- 
lutionary phases, that is. the populations of the most 
luminous stars which dominate the total integrated lu- 
minosity. Therefore once we explicit the relation be- 
tween the Wi coefficients and the most luminous stars, 
the values of luminosity ratios (i.e. colours) are also 
fixed. Since the proportionality relations between the 
stellar evolutionary phases that would be present in a 
cluster depend on the age of the cluster, in general 
Wi = Wi(t). This provides a way to estimate the age 
of the given cluster, for example by the comparison of 
different colours of the cluster with theoretical predic- 
tions. 

Note, however, that these relations do not allow us 
to obtain the total mass of the cluster unless the most 
luminous stars are also the ones that define the total 
mass (i.e. the unlikely case where the most massive 
stars are also the most numerous ones, a hypothesis 
that is ruled out by the observations). So what is the 
fraction of the total mass which is responsible for the 
total luminosity? Since we do not know the actual stel- 
lar mass distribution function of the cluster, we have 
to use a statistical method to describe how many stars 
with a given initial mass are expected in the cluster: 
the stellar initial mass function. For simplicity, we as- 
sume that all the stars in the cluster have been formed 
in a single star formation episode, and that there are no 
other episodes so that we do have a single stellar popu- 
lation (SSP). The integration of the stellar initial mass 



function over the mass range of initial masses that de- 
fines an evolutionary phase i, mi±dmi, provides the Wi 
coefficients that allow us to obtain the mass-luminosity 
ratio of stellar clusters as a function of age as 
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The inferred cluster initial mass is then obtained from 
this implicit mass-luminosity relation through a direct 
combination of the theoretical mean luminosity 



-^thco = y ] Wi £i , 
1=1 

with the inferred mass-luminosity ratio 



-^^infcrrcd ~ ^thco X [ -j- 



(7) 



(8) 



Strictly speaking, this direct comparison provides in 
fact the expected number of stars in the cluster < N >, 
and the expected mass < M > which corresponds to 
this expected number of stars. It is important to note 
that an incorrect age estimation also implies an incor- 
rect cluster mass estimation. 

Synthesis models provide this mass-luminosity re- 
lation for different ages and metallicitief0 for theo- 
retical clusters which contain an infinite number of 
stars. This deterministic method has been used, rather 
blindly, to clusters of any mass, and in this case the 
mass-luminosity relation is a simple function of the 
age and metallicity {M/L — f{t,Z)), and the funda- 
mental properties of clusters are inferred by statisti- 
cal tests such as fits. This result is also recovered 
as the mean value of the distribution of the possible 
mass-luminosities relations , (M/L = f (t, Z. M)), unde r 
a probabilistic framework ( Cerviiio fc Luridiana 20061 ). 
In this case, it is necessary to take into account the 
shape of the distribution and how it varies with the 
cluster mass when synthesis models results are applied 
to the analysis of real clusters. 



2 The distributed of the mass-luminosity ratio 
in stellar clusters and the initial cluster mass 
function 

The deterministic method is only valid, as seen above, 
in the limit of a very large number of stars populat- 



^Note that this metalhcity refers to the evolutionary tracks and 
not the metallicity in the stellar atmospheres, which may not be 
the same. 



Initial mass cluster distribution function 
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ing most, if not all, the evolutionary phases. The blind 
application of the method to small clusters results in 
wrong inferences have been made because the underly- 
ing assumptions are violated. A simple illustration is 
provided by Fig. [l] which shows the mass-luminosity 
ratio for different cases as a function of the cluster lu- 
minosity in the V band. The thick blue line at the top 
shows the evolution with age of the mass-luminosity ra- 
tio and the V magnitude for a 10'' Mg^ cluster obtained 
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SSP;10' M„ cluster 



from the SSP models provided bv lGirardi et al.l (|2002l ). 
The circle at the top of the line shows the position of a 4 
Ma-old cluster and the monotonic decrease in luminos- 
ity yields a monotonic increase in the mass-luminosity 
ratio, a property often used to infer fundamental prop- 
erties as described above. The thick black line at the 
bottom of the plot shows the position of 4 Ma-old indi- 
vi dual stars from the co rresponding isochrone provided 
by lGirardi et al.l (j2002l ). Note that, because of the use 
of the V band, the turn-off point appears to be the 
brightest with the smaller mass-luminosity ratio. The 
upper branch are the post-MS stars, while the lower 
branch provides the locus of the MS stars, down to 
very low luminosities and hence large mass-luminosity 
ratios. The shadow region in the middle of the plot is 
the result of 10^ Monte Carlo simulations for 4 Ma-old 
clusters using the same stellar initial mass function than 
the SSP models. The mass of each cluster is the result 
of a random sampling of a power-law initial cluster mass 
function with slope a = —1 covering the cluster mass 
range between 0.1 and 10^ Mq. The stellar initial mass 
function has been sampled randomly until the cluster 
mass of the cluster has been reached. 

Let us consider the mass-luminosity relation of in- 
dividual stars of a given age (bottom thick black line 
in Fig. [U isochrone of a 4 Ma-old population of sin- 
gle stars). Obviously, the possible values of the mass- 
luminosity depends on the particular properties of each 
star (its age, metallicity and mass; M/L = f{t, Z, M)). 
The figure shows that the range covered by the mass- 
luminosity relation of individual stars include all the 
possible mass-luminosity relations of SSP models. It 
also shows that the stellar mass-luminosity relation de- 
fines a natural limit of the mass-luminosity relations 
obtained in the Monte Carlo simulations of clusters. 

When stars are combined to describe stellar clus- 
ters (following a given stellar initial mass function) , the 
mass-luminosity relation gradually collapses to a single 
mass-luminosity ratio. The origin of this evolution in 
the distribution function of the mass-luminosity ratio 
is simply explained by the right-hand side of Eq. ([S]): 
the actual fraction of stars in a given evolutionary 
phase rii/Ni, does not coincide with the theoretical 
value Wi, but fluctuates around it following a multi- 
nomial distribution (see Cervino fc Luridiana ( 2006[ l 
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Fig. 1 Mass-luminosity relation in three sets of observ- 
ables : (1) Top thick blue line ; a 10^ Mq 4 Ma-old SSP 
model (circle) which slowly and monotonically evolves with 
age; (2) Bottom thick black line : single stars in a 4 Ma-old 
isochrone, the upper branch being the post-main-sequence 
stars, and the lower branch the main sequence; (3) Shaded 
area: 10*" simulations of a population of coeval 4 Ma-old 
clusters, whose initial mass function is a power-law with 
slope -1, with lower and upper mass limits of 0.1 and 10^ 
Mq respectively. Note at the bottom the sequences formed 
by binaries, triple and quadruple systems, etc. Clearly, the 
behaviour of the mass-luminosity ratio in this regime of 
small-mass clusters cannot be extrapolated from the evo- 
lution of the more massive cluster at the top. 



for a detailed discussion, andlCervifio fc Valls-Gabaud 



( 20031) :or a quasi-Poisson formalism). If Ni^Wi is the 



expected number of stars in the i-th evolutionary phase, 
the Ui value of real clusters will be distributed around 
it, producing variations with respect to the expected 
total luminosity of the cluster, but almost no vari- 
ation in its total mass. Equivalently, variations in 
the number of low-mass stars yield variations in the 
total mass, but not in the total luminosity. Obvi- 
ously, the dispersion in the mass-luminosity ratio will 
be larger for clusters which have a smaller number 



persions in n^l 
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20061). Note that the distributed 



nature of the mass-luminosity relation is a result of the 
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intimate composition of stellar populations of real stel- 
lar clusters. Its physical nature implies that it remains a 
distribution even in the case of perfect observations per- 
formed in perfect telescopes with perfect instruments 
with no statistical observational errors. 

Only when the number of stars in a cluster is large 
enough (i.e. the cluster is bright enough) the mass- 
luminosity ratio obtained by SSP models becomes a 
reliable, unique and well-behaved quantity. In other 
terms, the assumption of a mass-luminosity relation in- 
dependent of the cluster mass is only valid for massive 
clusters, typically with masses larger than 10'^ Mq. 

We want to stress that the main issue due to the in- 
complete initial stellar mass function sampling is that 
the proportionality between the actual evolutionary 
phases in the cluster at a given age ti, ni{ti)/nj{ti), 
differs from the assumed one in the synthesis models, 
Wi{ti)/wj{ti). Not only it may well be not fitted by 
the models, but it could also be close to the propor- 
tion Wi{t2)/wj{t2) that corresponds to a different age 
t2- For example, young clusters without massive stars 
(due to the sampling of the stellar initial mass function) 
are systematically best fit by models at older ages be- 
cause older clusters do not have massive stars. Under 
the usual assumption of full sampling, the sparse sam- 
pling of the IMF in these clusters is wrongly interpreted 
as a pure evolutionary effect. As the mass-luminosity 
ratio decreases with age, this effet translates into an 
overestimation of the initial cluster mass, producing a 
systematic bias in the cluster mass estimation. 

From another perspective, when sampling effects are 
present, there is more information on the properties of 
particular stars in the clusters (the effective temper- 
atures and luminosities of individual stars) but there 
is less information about the global properti es of the 
system (ag e and cluster in asses) . We refer to iBuzzoni 



obtain the most detailed information about the num- 
ber of stars in each evolutionary phase, that is, to 
analyse the colour-magnitude dia gram (i.e. the indr 
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for a more detailed analysis 
on the information that can be obtained from a stellar 
population through synthesis models. 

In the case of extreme sampling effects, the inte- 
grated light does not provide any information about 
the cluster, and accurate age or mass determinations 
can only be done taking into account the theoretical 
probability distribution functions that will produce a 
distribution of possible physical properties compatible 
with the observations. The range of physical properties 
will be larger when the number of stars in the clus- 
ter is smaller (the range of stellar mass-luminosity ra- 
tio is larger than the range predicted by SSP models), 
and implies an intrinsic loss of p recision in the global 
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vidual stars) of the clusters (e.g . iPellerin et al.l 12006 ; 
Hernandez fc Valls-Gabaudl 2008[ l. Unfortunately, the 
colour-magnitude analysis is not reliable for obtaining 
cluster masses, which are controled by low-luminosity 
stars. 

An alternative choice for a rough estimation of clus- 
ter ages is to look for signatures that are only present in 
a limited temporal range. In the case of young clusters, 
an example would be to look for young star signatures, 
such as Wolf-Rayet features or emission lines: the pres- 
ence of these signatures implies a young cluster, but the 
absence of these signatures does not imply an old clus- 
ter, just the absence of massive stars ! Again, this rough 
age estimation does not provide information about the 
cluster mass. 

2.1 A simple test for sampling effects identification 

The most trivial test to identify when sampling ef- 
fects are essential for an accurate analysis is to use the 
Lowest Luminosity L i mit ( LLL) method described in 
Cerviiio fc Luridianal (2004). The LLL implies that it 
has no meaning to compare a cluster with synthesis 
models (in a deterministic way) if the integrated lumi- 
nosity of the cluster is lower than the luminosity of the 
most luminous star included in the model. This simple 
statement restricts the deterministic use of synthesis 
models to young clusters with ma sses larger than a few 
10* M0 in the optical domain (|Cerviiio fc Luridiana 



20041) ■ which corresponds to a limiting magnitude of 
My = — 11 . In fact it i s just a common-se nse re- 
quirement: as an example, IZhang fc Faiil (Il999l ) reject 
point-like sources fainter than My = — 9 mag in the 
analysis of clusters in Antennae, since it is the luminos- 
ity of single luminous variable stars. However, based on 
the LLL requirement, not only "point-like" sources but 
all sources fainter than My = — 11 should be rejec ted in 
their analysis. For example. IPessev et al. ( 20081 ) show 
that young (200 Ma < t <1 Ga) clusters in the LMC do 
not fulfill the LLL requirements, and therefore the use 
of synthesis models within a deterministic framework is 
useless because it may yield wrong results. 

2.2 Implications for the initial mass cluster 
distribution estimation 



prope rties of the cluster (see 
20071 for a more extended discussion). The only 
way to estimate precise ages in this situation is to 



In a recent pedagogical paper. [Fall ( 2006[ l gives the re- 
lations between luminosity, mass and age distributions 
of young stellar clusters. A power-law luminosity func- 
tion for young clusters is directly related to a power- 
law initial cluster mass function, under the assumption 



Initial mass cluster distribution function 



5 



that the mass-luminosity ratio depends only of the age 
of the cluster. We have shown that this assumption is 
only valid for the case of massive clusters. The current 
controversy on the shape of the initial mass function 
of clusters, where small differences between a possible 
power-law or a log- normal distribution are important, 
depend crucially on clusters with masses around 10^ 
Mq. As we have shown, this mass range is below the 
limit of application of synthesis models in a determin- 
istic way, and a probabilistic framework is required for 
proper results, even though it implies an intrinsic loss 
of precision. 
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